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Introduction:  In  September  2007,  NRL,  in  part¬ 
nership  with  multiple  institutions,  undertook  a  com¬ 
bined  airborne  multi- sensor  remote  sensing  campaign 
and  in  situ  validation  effort.  The  experiment,  VCR’07, 
took  place  at  the  Virginia  Coast  Reserve  (VCR),  a 
National  Science  Foundation-funded  Long  Term  Eco¬ 
logical  Research  (LTER)  Site  on  the  Eastern  Shore  of 
Virginia.  The  study  area  comprised  an  1880  km2  region 
of  barrier  islands,  shallow  water  lagoons,  and  main¬ 
land  marsh  systems  (Fig.  4).  This  article  describes  the 
results  of  a  subset  of  experiments  conducted  during  the 
campaign  demonstrating  the  retrieval  of  soil  bearing 
strength  directly  from  hyperspectral  remote  sensing 
on  the  VCR  barrier  islands.  Bearing  strength,  or  “traf¬ 
ficability,”  is  a  key  parameter  needed  by  military  plan¬ 
ners  to  identify  littoral  penetration  points.  The  study 
also  developed  and  tested  new  methods  for  retrieval  of 
shallow  water  bathymetry,  another  important  param¬ 
eter  needed  by  amphibious  craft  during  landing. 

Background:  Hyperspectral  imaging  (HSI)  sensors 
record  the  reflected  solar  radiation  from  land  and  water 
and  have  been  used  to  retrieve  important  information 
in  a  variety  of  applications  including  precision  land- 
cover  mapping;  in-water  retrievals  such  as  bathymetry, 
bottom  type,  and  suspended  constituents;  retrieval  of 
biophysical  and  geophysical  parameters  on  land;  and 
detection  and  determination  of  man-made  structures 
and  objects.  HSI  sensors  are  unique  in  having  a  large 
number  of  narrow,  contiguous  spectral  channels.  They 
have  sufficient  spectral  resolution  to  identify  different 
surface  characteristics  based  solely  on  spectral  signa¬ 
tures,  and  allow  mapping  of  retrieved  quantities  with  a 
high  degree  of  precision. 


Airborne  Data  Acquisition:  For  VCR’07,  NRL 
mounted  three  remote  sensing  instruments  in  a 
de  Havilland  Twin  Otter  aircraft  (Fig.  4):  a  CASI- 
1500  visible  near-infrared  (VNIR)  hyperspectral 
camera  operating  in  the  0.38-1.04  micron  spectral 
range,  a  Surface  Optics  hyperspectral  short-wave 
infrared  (SWIR)  camera  operating  in  the  0.9- 1.7 
micron  range,  and  a  single- channel  mid- wave  infrared 
(MWIR)  camera  operating  in  the  3-5  micron  range. 
This  article  focuses  on  the  results  obtained  using  the 
CASI- 1500. 

Bearing  Strength  from  Hyperspectral  Imagery: 

Soil  bearing  strength  depends  on  a  number  of  key 
factors  including  soil  composition,  grain  size,  and 
water  content.  HSI  sensors  can  discern  these  prop¬ 
erties  for  the  surface  layer  of  a  substrate:  HSI  soil 
composition  mapping  has  been  well  documented 
in  the  literature;1  some  studies  have  demonstrated 
the  ability  to  model  grain  size  from  HSI;2  and  liquid 
water  absorption  features  visible  at  several  parts  of  the 
reflectance  spectrum  allow  quantitative  estimates  of 
moisture  content. 

To  determine  the  relationship  between  HSI  surface 
data  and  the  potential  bearing  strength  of  the  substrate 
overall,  the  VCR’07  team  undertook  in  situ  spectral 
and  geotechnical  measurements  to  characterize  bearing 
strength,  moisture  content,  grain  size  profile,  and 
corresponding  spectral  reflectance  of  different  typical 
substrates  (Fig.  5).  These  measurements  were  used  to 
develop  a  spectral  look-up  table  for  bearing  strength, 
which  was  then  used  to  produce  maps  of  substrate 
bearing  strength  directly  from  the  CASI  hyperspectral 
imagery.  Figures  6  and  7  show  examples  of  in  situ  data 
retrieval.  Figure  8  shows  the  CASI  retrieval  product. 
The  validation  efforts  showed  that  although  HSI  sees 
only  the  surface  reflectance,  this  layer  provides  diag¬ 
nostic  information  to  indicate  an  estimated  bearing 
strength  of  the  substrate. 

Shallow  Water  Bathymetry:  At  visible  and  near- 
infrared  wavelengths,  reflectance  from  the  water 
column  is  determined  by  a  variety  of  factors  including 
water  depth,  bottom  type,  and  the  presence  of  sus¬ 
pended  constituents  such  as  color-dissolved  organic 
matter  (CDOM),  suspended  sediments,  chlorophyll, 
and  phytoplankton.3  Thus  the  general  problem  of 
retrieving  depth  as  well  as  other  water  properties 
has  been  approached  using  spectral  look-up  tables4 
in  which  a  forward  radiative  transfer  model  such  as 
Hydrolight5  is  executed  repeatedly  with  varying  water 
column  properties,  depth,  and  bottom  types.  To  be 
comprehensive,  these  look-up  tables  must  be  large  and 
may  need  to  be  tuned  to  specific  coastal  types  because 
bottom  types  and  water  properties  may  vary  signifi¬ 
cantly  with  coastal  type. 
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Hog  Island 
Wreck  Island 
Smith  Island 
Indiantown  Marsh 


FIGURE  4 

Left:  Twin  Otter  aircraft  and  the  three  NRL  sensors  (CASI-1  500,  Surface  Optics  SWIR,  and  mid-wave 
IR)  during  VCR'07.  Right:  Composite  of  VCR'07  CASI-1 500  quicklook  images,  with  arrows  indicat¬ 
ing  locations  of  primary  bearing  strength  retrieval  validation  sites. 


In  very  shallow  water  (<  2  m  depth),  there  is  a 
range  of  wavelengths  in  the  near- infrared  for  which  the 
dominant  factors  determining  reflectance  are  bottom 
type  and  depth,  with  water  column  constituents 
playing  a  secondary  role.  During  VCR’07,  we  measured 
reflectance  as  a  function  of  depth  for  various  typical 
bottom  types,  and  determined  the  optimal  wavelengths 
for  which  a  simpler  model  involving  a  library  of  regres¬ 
sions  for  each  bottom  type  could  optimally  be  applied 
in  the  very  shallow  limit.  One  of  the  optimal  wave¬ 
lengths  for  this  occurs  near  810  nm,  a  local  minimum 
in  liquid  water  absorption6  and  thus  a  local  maximum 
in  reflectance  (Fig.  9).  The  results  shown  in  Fig.  10 
demonstrate  the  retrieval  of  very  shallow  depth  in  the 
vicinity  of  one  of  the  VCR  barrier  islands  using  this 
approach.  Kinematic  GPS  data  taken  in  the  intertidal 
zone  were  used  to  demonstrate  the  high  accuracy  of  the 
retrieved  very  shallow  bathymetry  product. 

Conclusions:  VCR’07  successfully  demonstrated 
new  retrievals  from  hyperspectral  imagery.  These 
included  a  novel  bearing  strength  map  suitable  for 
trafficability  analysis  in  a  barrier  island  coast  type.  A 
simplified  approach  to  retrieving  bathymetry  from  HSI 
in  the  very  shallow  limit  was  also  validated.  Opera¬ 
tional  lidar  systems  such  as  Scanning  Hydrographic 
Operational  Airborne  Lidar  Survey  (SHOALS)  do 
not  produce  reliable  depth  retrievals  in  very  shallow 


waters.  In  joint  HSI/lidar  platforms  such  as  the 
Compact  Hydrographic  Airborne  Rapid  Total  Survey 
(CHARTS),  the  new  approach,  therefore,  could  fill  a 

gap  in  bathymetric  retrieval. 

[Sponsored  by  ONR  and  the  National  Geospatial- 

Intelligence  Agency] 
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FIGURE  5 

In  situ  calibration/validation  instruments:  Analytical  Spectral  Devices  (ASD)  full-range  spectrometer;  lightweight 
deflectometer  (LWD),  which  measures  dynamic  deflection  modulus  (bearing  strength);  dynamic  cone  penetrometer 
(shear  strength);  and  soil  grab  sampler  (for  laboratory  grain  size  profile  and  moisture  analyses). 


07HIA-1 ,  Evd  =  24.3  MN/m2,  Moisture  =  20.36% 


0.10 
8  0.05 

c 
as 

I  0.00 

DC 

-0.05 

-0.10 

500  1000  1500  2000  2500 

Wavelength  (nm) 


CD 

o 

c 

03 

o 

0) 

Q 


DC 


07HIA-13,  Evd  =  6.3  MN/m2,  Moisture  =  1.03% 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

500  1000  1500  2000  2500 

Wavelength  (nm) 


FIGURE  6 

In  situ  water  content /spectral  measurements.  ASD  reflectance  spectra  for  (left)  sand  near  the  shore¬ 
line  and  (right)  sand  near  the  arid  backdune.  The  shoreline  site  with  a  moderate  measured  moisture 
level  shows  liquid  water  absorption  features  in  the  spectrum;  the  backdune  site  with  a  low  measured 
moisture  level  has  a  spectrum  absent  of  these  liquid  water  absorption  features.  The  corresponding 
bearing  strength  measured  by  the  LWD  was  high  at  the  shoreline  site  and  low  at  the  backdune  site. 
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Wreck  Island  Transect  C,  Position  8  (WIC-8) 
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quartz  sand, 
no  silt. 


FIGURE  7 

Correlation  of  in  situ  grain  size,  moisture,  and  bearing  strength  measurements.  Scatterplot  of  percent  moisture  vs 
LWD-measured  dynamic  deflection  modulus  (bearing  strength)  for  beach  shoreline  and  backdune  samples.  Pres¬ 
ence  or  absence  of  silt  leads  to  dramatically  different  bearing  strength  when  moderate  moisture  is  present.  Since 
silt  is  a  significant  contributing  factor  to  bearing  strength,  we  looked  for  diagnostic  features  that  might  indicate  the 
presence  or  absence  of  silt  quantitatively.  Correlation  of  grain  size  and  ASD  spectral  data  (not  shown)  provided 
evidence  for  the  presence  of  diagnostic  VNIR  spectral  features  for  grain  size. 


FIGURE  8 

Demonstration  of  a  bearing  strength  retrieval  product.  Left  to  right:  CASI-1500  HSI  scene,  Hog  Island,  VA;  clos¬ 
est  matching  spectrum  in  the  spectral-geotechnical  look-up  table  (LUT);  and  resulting  retrieved  bearing  strength 
estimate  in  MN/m2. 


226 


2009  NRL  REVIEW 


Reflectance 


REMOTE  SENSING 


—  water_90cm 

—  water_86cm 
water_83cm 
water_74cm 
water_70cm 

—  water_60cm 

—  water_56cm 

—  water_45.5cm 
water_37.5cm 
water_29cm 
water_21cm 
water_1 3cm 
water_8cm 

—  muddy_bottom_ 
tidal_creek_bank 


(D 

CL 


FIGURE  9 

Top:  In  situ  spectrometry  site  on  salt  marsh  tidal  creek  bank  with  muddy  bottom,  at  high  and  low  tide,  Parramore  Island,  VA.  Bottom 
left:  In  situ  spectral  reflectance  vs  depth  profiles.  Bottom  right:  Regression  at  the  810  nm  feature. 


2009  NRL  REVIEW 


227 


REMOTE  SENSING 


FIGURE  10 

Shallow  water  bathymetry 
product.  Left:  NRL  CASI-1500 
image  of  Wreck  Island,  VA. 
Right:  Retrieved  depth  using 
the  810  nm  regressions  for 
each  bottom  type  determined  in 
preprocessing  from  the  spectral 
libraries;  depth  is  quantized 
in  the  color  scale,  but  actual 
retrieval  is  continuous. 
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